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O . Abstract 



XEFT is a low energy effective theory for the X(3872) that can be used to systematicahy analyze 
the decay and production of the X(3872) meson, assuming that it is a weakly bound state of 



{2L|' charmed mesons. In a previous paper, we calculated the decays of X(3872) into Xcj plus pions using 

. a two-step procedure in which Heavy Hadron Chiral Perturbation Theory (HH^PT) amplitudes 

\ are matched onto XEFT operators and then X(3872) decay rates are then calculated using these 

XjT^ ' operators. The procedure leads to IR divergences in the three-body decay X(3872) — )■ XcJT^t^ when 

CN ■ virtual D mesons can go on-shell in tree level HHxPT diagrams. In previous work, we regulated 

o 

these IR divergences with the D*^ width. In this work, we carefully analyze X(3872) — )■ XcJT^^ 
and X(3872) — t- XcJT^t^ using the operator product expansion (OPE) in XEFT. Forward scattering 
amplitudes in HHxPT are matched onto local operators in XEFT, the imaginary parts of which are 
^ . responsible for the decay of the X(3872). Here we show that the IR divergences are regulated by 

I the binding momentum of the X(3872) rather than the width of the D*^ meson. In the OPE, these 

IR divergences cancel in the calculation of the matching coefficients so the correct predictions for 
the X(3872) — )■ Xcitttt do not receive enhancements due to the width of the D*^ . We give updated 
predictions for the decay X(3872) — )• Xcitttt at leading order in XEFT. 
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''' Electronic address: melien@pliy.duke.edu 
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I. INTRODUCTION 

The X(3872) [l-3| is most likely a molecular bound state of neutral and D*^ mesons 
(and their antiparticles) . Recently we developed an effective field theory for the X(3872) 
called XEFT which can be used to systematically analyze its properties under this assump- 
tion j4|. This theory contains nonrelativistic vr^, D", and D*^ mesons as its degrees of free- 
dom. It has been used to rederive effective range theory predictions for X(3872) — t- D^D^tt^ 
and to calculate the leading corrections to this process [4|, to calculate elastic scattering of 
charm mesons off the X(3872) M, the cross section for 7r+X(3872) d*+D*^ 6 and to 
calculate decays of X(3872) to quarkonia [7|, 18|. Decays to final states with Xcj are interest- 
ing because heavy quark symmetry predicts the relative rates for decays to different Xcj and 
therefore the relative rates can be used to test the molecular interpretation of the state {q]. 

The purpose of this paper is to carefully revisit the calculation of three-body decays, 
X(3872) — 7- XcjTTvr, first performed in Ref. |7|. In the case of decays to Xci, Ref. {3] found 
double poles coming from static D meson propagators that can go on-shell in certain regions 
of three-body phase space. These double poles lead to divergent phase space integrals in the 
total decay rate. Since the divergence is associated with an internal D meson going on-shell 
we will refer to this as an IR divergence. Ref. ^ regulated these double poles by performing 
a matching calculation with D*^ having a complex nonrelativistic energy A + iTo*" /2, where 
A is the hyperfine splitting between the and D*'^ and r/)*o ^ 68keV is the width of the 
D*^. When the double pole is regulated this way the decay rate is enhanced by a factor 
D*o ~ 10^ where E^^ is the typical energy of the pion in the decay. Therefore, 
predictions for decay rates for X(3872) — )■ Xcitt^tt^ in Ref. \l\ were found to be quite large. 
In this paper, we argue that the IR divergence is not regulated by the width of the D*^ 
but by larger mass scales. We show how to determine how to properly handle these IR 
divergences in these calculations and give the correct prediction for the three-body decay 
X(3872) —7- XciTTTT. The issues addressed in this paper should be relevant to other decays of 
the X(3872) as well as the decays of other states that could be molecular bound states of 
hadrons, e.g., the recently discovered states ll|]. (For an XEFT-like treatment of these 
states, see Ref. jl2|)- 

XEFT is supposed to be valid near (within a few MeV) of the D^D*^ threshold and 
other nearby thresholds, such as the D^D*^ threshold which is about 8 MeV away are 
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not explicitly included in the calculations mentioned above. A recent calculation of the line 
shape in the vicinity of the X(3872) in Ref. [3] shows that this is a very good approximation. 
It would be interesting to revisit many XEFT calculations with charged mesons included 
explicitly in the theory but we will not do so in this paper. 

In order to understand the nature of the IR divergences we consider an alternative ap- 
proach to calculating the decay rate X(3872) — )■ Xcjttvt. We begin by reviewing how the 
hadronic decays to quarkonia plus pions were calculated in Ref. |7|, using X(3872) — )■ Xcjtt^ 
to illustrate the technique. The underlying short-distance process is D^D*^ — t- Xc.j'^^-^ By 
short- distance we mean that the process occurs on length scales much shorter than the typ- 
ical separation of the constituents of the X(3872). This length scale is not well known since 
the binding energy of the X(3872) is not accurately determined, but is expected to be large 
compared to typical hadronic scales because the X(3872) is very weakly bound. The wave 
function of the D^-D*^ at long-distances is 

g-7'' 

i^DD-ir) , (1) 



where 7 = y/2n£)£)*B, where Udd* is the reduced mass of the charmed mesons and B is the 
binding energy. From the currently measured masses of the X(3872),D°, and D*^ [13] we 
know that B = 0.42 ± 0.39 MeV. For this range of binding energies, the root mean squared 
separation is rx = he/ {\/2'y) = 4.9l}^4^ This is much larger than the size of the 
and D*^ meson so the charm and anticharm quarks in the X(3872) are well separated most 
of the time. To coalesce into a quarkonium the D mesons must come together to within 1 



fm, a sca^ 
theorem 



e much shorter than hc/'j. One then expects to be able to derive a factorization 



14j-ll6l| which schematically takes the form 

r[X(3872) ^ L.H.] oc {tpoD* (0) | V[D°D*° + c.c ^ L.H.] , (2) 

similar to what is found for the decay of quarkonia and other nonrelativistic bound states. 
However, the wavefunction in Eq. ([T]) diverges as r — i- so Eq. ([2]) is not quite right. The 
description of the X(3872) as a two-body bound state surely fails as r — )■ 0. New channels 
such as J/ipp and D^D*^ + c.c. open up and at some point a hadronic description in terms 
of mesons ceases to work entirely. In XEFT this is dealt with by requiring that the cutoff. 



^ Here and below charge conjugated processes are implied. 



3 




FIG. 1: Feynman diagram contributing to the X(3872) — )■ XcJT^^ decay amplitude. The circle with 
cross represents the interpolating field for the X(3872). 

A, be taken not too much larger than a few times 7 so these effects are integrated out and 
included as short distance operators in XEFT. Thus we expect \ipDD*{0)\'^ to be replaced by 
the appropriate short-distance matrix element in XEFT. 

The factorization theorem of Ref. j?! is derived as follows: the amplitude for the transition 



D D — )■ XcjTT is computed in Heavy Hadron Chiral Perturbation Theory (HHxPT) 



19| . The amplitude contains contributions from contact interactions and tree level graphs 
with virtual D mesons. For the two-body decays, the D mesons are off-shell by E^^ + cA, 
where c = or ±1 depending on the channel. Since E^^ ranges from 305 to 432 MeV in 
these decays and A = 142 MeV, the virtuality of the D meson is large compared to the 
relevant scales in XEFT and the virtual D mesons should be integrated out. The amplitude 
for D^D*^ — 7- Xcqt^^, for example, is reproduced in XEFT by a local operator coupling the 
four fields: 

Here are the fields for D*^ meson fields and P are the fields for the D mesons. The 

function oI-^tt) is calculated in Ref. and contains the energy dependence of the virtual 

D mesons. The decay of the X(3872) is then calculated using the Lagrangian in Eq. Q. 

The relevant Feynman diagram, depicted in Fig. [H gives the result (for arbitrary J) 

7(A pDs - 7)^ m^^j Vl 
27r mx 727r/2 

where ao = 4, ai = 3, and 0^2 = 5, the functions C^^j{Et,) are given in Ref. [7|, and the PDS 
subtraction scheme is used to regulate the divergent integral appearing in the calculation of 
the rate. The first term in Eq. (jl]) is equal to the following matrix element in XEFT 

7(APDS^- 1? ^ 1 ^ |(o|-Le,(A) (V^ P + V^ P)|X(3872, A))p . (5) 

A 



r|X(3872) ^ x.jA = iiliS^^ ^ a, \CUE.f . (4) 



Thus we find that the decay rate factorizes directly into a matrix element which plays the role 
of the wavefunction at the origin squared times the remaining factors which are calculable 
in HH^PT. The matrix element is not calculable in XEFT, and the sensitivity to short 
distance physics is evident from the quadratic dependence on the cutoff. 

In next section of the paper, we briefly discuss power counting in HH^PT and XEFT. 
In the following section, we will show how to rederive the previous results for two-body 
X(3872) decays using the operator product expansion (OPE). In the next section of the 
paper, we will apply the OPE to the three-body decays of the X(3872) where we reproduce 
our previous results for the decays X(3872) — )■ Xcj^rvr for J = and 2, and provide the 
correct treatment of IR divergences for the case J = 1. 



II. POWER-COUNTING 



Before applying the OPE in matching HHxPT onto XEFT, we explain the power counting 
in XEFT and HHxPT. In HHxPT the expansion parameter is Q/A^, where = 47r/^ ~ 
1 GeV. Here Q is either rriT^ or p or E, where p (E) are either a loop or external momentum 
(energy) that are taken to be of order m^. In XEFT the scales m^. A, and similar scales have 
been removed from the theory and reside in short-distance coefficients. As a result the short- 
distance coefficients scale with powers of m.,^/Ay (where we have chosen to replace Q with 
for bookkeeping purposes). In Ref. j4], XEFT is formulated as an expansion in Qx/Ax 
where Qx ^ Pw ^ Pd ^ Pd 1 = a/A^ — = 44 MeV, ^ and the unknown scale. Ax, 
sets the range of the convergence of XEFT and is expected to be of order m^. The XEFT 
operators can have a complicated energy or momentum scaling since both A^ and scales 
of order rrij^ can appear in the denominator of XEFT operator coefficients. However, the 
scaling of XEFT operators becomes simple when all momentum scales Qx are considered as 
being mj^VT,, where v-,, ^ 1/3, for the purposes of power counting. The pion kinetic energy, as 
well as energies in loop integrals, are OlQ'j^/mj^), and therefore are counted as m^rf^. In this 
formulation of the power counting, powers of from momenta and energy in the numerator 
cancel powers of in the denominator so that a generic XEFT operator scales as powers 
of (m7r/A^)™'f where the integers m and n positive. In this formulation, calculations in 



The scale n appears in potential pion exchange graphs in XEFT 
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FIG. 2: The operator product expansion for D^D*^ — )• Xcj^^^ — ^ D^D*^. The forward scatter- 
ing matrix element in HHxPT on the left-hand side is matched onto products of short-distance 
coefficients, Aj, and XEFT operators on the right-hand side. 

XEFT are a simultaneous expansion in 111^^/ and v-,^. Examples will be given below. 



III. TWO-BODY DECAYS AND THE OPERATOR PRODUCT EXPANSION 

Strictly speaking, matching amplitudes in HHxPT onto local operators in XEFT as in 
Eq. ([3]) does not adhere to the EFT paradigm because in two-body decays the energies of 
the Xcj and 7r° in the final state are too large for the power counting of XEFT. In the decay 
X(3872) ^ Xcjvr°, the energy of the pion is 433 MeV, 347 MeV, and 305 MeV, for J = 0, 1 
and 2, respectively. The pion is relativistic and therefore is not properly treated in XEFT. 
Likewise the momentum of the Xcj exceeds the size mandated by XEFT power counting. 
These particles should really be integrated out of XEFT and their effects incorporated into 
local XEFT operators. 

The rigorous approach to X(3872) — )■ Xcjt^^ is to use the OPE. In this approach the 
forward scattering matrix element for D^D*^ — t- D^D*^ via an intermediate Xcjt^^ state 
is calculated in HHxPT, and matched onto products of short-distance coefficients and lo- 
cal operators in XEFT. This is illustrated in Fig. [2l After matching, the decay rate for 
X(3872) — 7- XcjTT*^ can be calculated in XEFT. This approach is conceptually similar to 
how short- distance annihilation decays of quarkonia are treated in Non-Relativistic QCD 
(NRQCD) 2o|. In that theory, full QCD is used to calculate the matrix element squared 
for the processes cc — > gg,ggg, for example. The matrix elements squared in full QCD 
determine the imaginary parts of coefficients of four-quark operators in NRQCD and the 
imaginary parts are responsible for the decay of the quarkonium state in NRQCD. 

The HHxPT forward scattering matrix elements for D^D*^ — t- D^D*^ via intermediate 
Xcj7i^° states are shown in Fig. |3]for the Xco (top row), Xci (middle row), and Xc2 (bottom 
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FIG. 3: Diagrams for the forward scattering amplitude for D°D*° XcJT^^ D^D*^ in HHxPT. 
The top, middle, and bottom rows show the diagrams involving the XcOi Xci-, S'ld Xc2-, respectively. 

row). The amplitudes for these diagrams are 



- Kf ) / 



^ '' d'^i — |/p 1 rriy, 1 mp*mp 



(27r)4 (£0)2 e-ml + ie mx Qi - £° + ie 2 



9 V A / i (27r)4 + A)2 P-ml + ie mx - ^° + ie 2 
where Qj = mx — f^xd- The factor of Tn^^j/mx is associated with each Xcj propagator. ^ 
To determine the two-body decay rate, we will need the discontinuity of the expressions 
above. The discontinuity can be found by cutting the both the pion and XcJ propagators, 
which amounts to the following replacement rules: 

^ ^ -in 5{VLj - e^) . (8) 



- £0 + ie 



^ This factor can be understood by starting with a fully relativistic kinetic term for the XcJ fields, 

d^^xljd^XcJ - ml^jXljXcJ , 

then making the field redefinition 



The factor of ^Im^^ j ensures proper normalization of the non-relativistic Xc.i field and the factor mx ~ 
todo + m£)»o in the exponent appears because all energies are measured relative to the D'^D*'-^ threshold. 
The factor m-^^j/mx could consistently be set to 1 within the accuracy of our calculation, but we keep it 
to obtain expressions for decay rates that are identical to our earlier calculations. 
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Using these rules we find: 



2z ^ 8n3\Uj Elmx 2 



where in each expression E^, = Vtj and p^, = ^^Q'j — m^. 

Next we need the local XEFT operators which we will use in the matching: 

5£ W = -J2 Af\v'P + V'PyiV'P + V'P) , (10) 
J 

where the subscsript n refers to the number of particles in the final state and in this section 
we are focusing on n = 2. The amplitudes for D^D*^ D^D*^ for these operators are 

-i^Af AmD*mr), (11) 
J 

and matching onto the HH^PT amplitudes in Eq. (|6]) gives 



4(2) _ 
J\q — 






A? = 






A'? = 







(12) 

The matching coefficients scale as 

/I (2) 2"^^ /1Q\ 

X 

in XEFT according to the power counting discussed earlier, since i?^ ~ ~ A ~ O(m^). 
Note that because the HHxPT forward scattering amplitudes match directly onto local four 
fermion operators in XEFT there is no suppression due to factors of ^v• 

The XEFT decay rate for X(3872) coming from final states with n particles {n — 1 pions) 
is given by 

2ImS („)(-Ex) 
~ ReS'(-Ex) ' 
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FIG. 4: The leading contribution to S (2)(— -Ex)- 

where Tj{—Ex) is the sum of Cq irreducible graphs contributing to the two-point function of 
the X(3872) interpolating field, evaluated at the X(3872) pole. S -Ex) is the contri- 
bution to S(— -Ex) from Cq irreducible graphs with one insertion of v4j"\ The leading order 



contribution to ReS'(— -Ex) was calculated in Ref. j4|: 

ReE'(-E;x)LO = ^. (15) 
zn-y 

The leading contribution to S (2)(— -Ex) is given by the diagram in Fig. HJ which evaluates 

to 

^(2) 

E^(2)(-i?x)LO = -7^. (16) 
Using these results in Eq. (fl^ . the leading order XEFT decay rate for X(3872) — )• XcJ^r^ is 

F^O[X(3872) -> XcjA = -^^ImAf , (17) 

where the extra factor of two comes from including not only the (V^PyV^P operator, but 
also the {V'PyV'P, {V'PyV'P, and {V'PyV'P operators, and a factor of 1/2 from the 
squaring the factor of 1/a/2 in the X(3872) wavefunction. This expression makes it clear 
that we are only interested in ImAj = DiscAj/(2i), which at leading order can be read off 
of Eqs. ([9]) and ( IT2l) . Inserting these into Eq. (ITTl) gives 

F,,,,[X(3872) ^ Xcivr ] - ^^^^^^^^ 

F,,,,[X(3872) ^ xc.n ] - 7^ (:^;;^ • (18) 

This agrees with the results of our previous paper 
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FIG. 5: HHxPT forward scattering matrix elements for D^D*^ — )• D^D*^ via intermediate XcJT^^^^^ 
states, with J = 0, 2. The dark solid line is the x state. 

IV. THREE BODY DECAYS AND THE OPERATOR PRODUCT EXPANSION 

In this section we apply the OPE to the three-body decays of the X(3872), which is 
straightforward for final states with a Xco or Xc2- The calculation of these decays is completely 
analogous to the two-body decays of the previous section. The decay rate for the three- 
body decays under consideration is given by Eq. (|T4l) . The matching coefficient Aj needed 
to calculate ImS (3)(— -Ex) is determined by matching HHxPT forward scattering matrix 
elements for D^D*^ — t- D^D*^ via intermediate Xcj'^^^c^ states onto XEFT . The diagrams 
involving Xco or Xc2 are shown in Fig. |5l We calculate the imaginary parts of these graphs 
using the optical theorem to determine Im Aj . The decay rate is then computed in terms 
of Im A J by evaluating the diagram in Fig. HI The final results for the LO decay rates for 
X(3872) ^ Xcjvr°vr°, J = and 2, are 



r[X(3872)^XcovrVi 

•^Qo—rriTr 

dEi 



27r7 



/i|,^.Co(ApDs)2 27 mx 



dE2 S{Ei + E2- no)plpl F'\Ei, E2 



(19) 



where 



F(Ei, E2) = {E2 — El 



1 



iEi + A){E2 + A){Ei + E2 + A) 

3 A 



+ 



El E2 {El + E2-A) El E2 {El + A) {E2 + A) 



(20) 
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and 



r|A-(3872) = , ^7, '-^ (21) 

(•02— m,^ pQ,2~mT, 

X dEi dE2 5{Ei + E2- ^2) 



CiEi, E2f + C(Ei, E2) E2) - D{E2, E,)) 

+D{Eu E2f + D{E2, E^f + E2) D{E2, E,) 



where 



C{E,,E2) 



E2 — El 



+ A)(E2 + A)(Ei + E2 + A) 



Here, Ei{pi) refers to the energy (three- momentum) of one of the 7r° and Pi = \pi\. The 
partial decay rates are symmetric under 1^2. In these decays A is equal to the neutral 



hyperfine splitting, A = m£)*o — m^io 







42 MeV. These expressions for the decay rates 



are identical to those calculated in Ref. [7|, except for how the phase space integration is 
performed. Ref. Q] calculated the rates with the exact relativistic phase space for the Xcj 
and pions. In the present paper we drop terms of order uIt^/Mx in the argument of the 
energy-momentum conserving delta-function in three-body phase space. ^ This considerably 
simplifies the calculation of three-body phase space, which is now a one-dimensional integral 
due to the delta-functions in Eqs. ( fT9l) and (12T1) . The rates calculated with this approxima- 
tion for the phase space agree with the rates in our previous calculation to an accuracy of 
10-15%. 

Next we turn our attention to the problematic decay rate for X(3872) — )■ Xci7r°vr°. The 
HHxPT diagrams for the forward scattering amplitude for D^D*^ — )■ D^D*^ with an inter- 
mediate Xci^r'^vr^ state are shown in Fig. [6l One would like to directly match these diagrams 
onto a local XEFT operator and obtain the coefficient Im A\ to calculate the decay rate 
as in our previous calculations. However, there are diagrams in which the internal D meson 
propagators can go on-shell. Both internal D meson propagators in Figs. [6^), b), e), and f) 



^ Ref. states that using the fuUy relativistic phase space is necessary because expanding in m^/Mx in 
the argument of the energy /momentum conserving delta- function would lead to unconstrained momentum 
integrals, but this is incorrect. 
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FIG. 6: HH^PT forward scattering matrix elements for D^D*^ 
states. The dark solid line is the Xci- 



via intermediate Xci^r vr 



0^0 



can go on-shell, leading to double poles in the three-body phase space integration. These 
double poles require a regulator to obtain finite answer for the total decay rate. The internal 
D meson propagators on the left-hand sides of Figs. Eb), d), g), and h), and on the right- 
hand sides of Figs. [6]), k), o), and p) can also go on-shell, which leads to single poles in the 
three-body phase space. These single poles can be dealt with a principal value prescription. 
The imaginary part of the diagrams in Fig. |6] is given by 



n*o 



(23) 



where 



iO n*0 



3 



(24) 



mx 



dEi 



dE2 6{Ei + E2-ni)pipi 



?,A{E,,E2f-2A{E,,E2) 



+ 2 



El E2 El E2 {El + E2 



+ 



El + El 



ElEl{Ei + E2Y 



+-A{Ei,E2 



Bi{Ei,E2) ^ Bi{E2,Ei] 



2 

+3 



El- A E2- A 
Bi{Ei, E2) Bi{E2, El) ^ Bi{Ei,E2) 1 
^1 - A E2-A 
Bi{Ei,E2) 1 1 



Bi{E2,Ei) 1 



Ei-A E2E1 + E2 
Bi{Ei,E2Y Bi{E2,EiY 



El — A El El + E2 
Bi{E2,Ei) 1 1 
E2-A E1E1 + E2 



E2 — A E2 El + E2 



{Ei-A)^ {E2-AY 
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FIG. 7: Long distance XEFT contributions to the X(3872) — t- Xci'^^'^^ decay rate coming from 
diagrams with two on-shell D meson propagators (on the left) or only one on-shell D meson 
propagator (on the right). 



where 



El E2 {El + E2) {El + A) {E2 + A) ' 

«'(^"^^) = ^ + irT^- 

The terms in the curly brackets are increasingly singular when the D meson goes on-shell, 
which occurs for either Ei ^ A or E2 ^ A. Note that it is not kinematically possible to 
have both Ei and E2 go to A at the same time. In Ref. [3], the integrals were rendered 
finite by replacing {Ei - Af {Ei - Af + V^/A, where T is the width of the D*°. This 
yields a finite result that diverges in the F — )■ limit. In what follows we will make a similar 
replacement to render the integral in Eq. finite, but will interpret F as a regulator rather 
than the physical width of the D*^. 

Because there are contributions to the graphs in Fig. [6] in which two D mesons and the 
pion can go on-shell, these graphs cannot be directly matched onto a local XEFT operator. 
There are also long distance XEFT contributions which must be taken into account; these 
are shown in Fig. [71 Only after subtracting these long distance contributions from Eq. ( 12^ 
can the remainder term be matched onto the product of a local XEFT operator and the 
short distance coefficient A^i^ in Eq. ( ITOj) . The diagrammatic equation for the matching 
coefficient A^^-* is shown in Fig. |8l where the blob on the left represents all of the HHxPT 
forward scattering amplitudes in Fig. |6l The contact interaction represented by the solid 
square in Fig. UK) can be obtained by matching the HHxPT forward scattering matrix 
element for DD — t- DD via an intermediate Xci^ state onto a local XEFT operator of the 
form 

B{DD)\DD) . (26) 
13 



FIG. 8: The matching calculation for the coefficient A\ . The blob on the left represents the sum 
of all diagrams in Fig. [6j 

We will determine Im B below. There is also another contact interaction in Fig. [7b) which 
is obtained by matching the HH^PT diagrams for D^D^* — ?■ D^D^-k^ . We will not need the 
explicit form of this operator. 

According to the XEFT power counting discussed in Section II, the diagram in Fig. [7^) 
scales as ((y'^m^f,r)/A^, which is f,r ~ 1/3 suppressed relative to ImAj"'' (n=2,3). This 
diagram has two D meson propagators that can go on-shell, so it is IR divergent and therefore 
must be included in the matching calculation. We will see below that the IR divergence in 
Fig. [7|i) will cancel the IR divergence in Eq. ( 12^ . rendering the coefficient ImA]^ finite. 
The power counting gives the correct estimate for the size of the remainder of the diagram, 
however, this remainder does not contribute to the X(3872) Xci7r°vr° decay width. The 
diagram in Fig. [7b) only has a single on-shell D meson propagator and scales as {glrn^v^) / A^. 
Thus this contribution is suppressed by Vt, relative to the IR finite part of Fig. [7^) and by 

relative to the coefficient A^^\ Furthermore, it is IR finite, i.e. it does not diverge as 
r — )■ 0, so we will drop this graph in the matching calculation as well as the calculation of 
the decay rate. 

Next we will calculate the loop diagram in Fig. [7^). We first calculate \m.B which is 
determined by the matching diagrams shown in Fig. [9] The HHxPT forward scattering 
matrix element for DD — )■ DD has the form 

= A\SS,^r^ml (27) 

so 

B = -lA^^,^r'- (28) 
The explicit expression for the HHxPT forward scattering matrix element calculated from 
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FIG. 9: Matching the HHxPT forward scattering matrix element for DD — )• DD via an interme- 
diate XciTi" state onto a local XEFT operator. 



the diagram on the left-hand side of Fig. [9] is 



,AgYi f dH |£f 1 1 



""^ j (277)4 (£0 _ A)2 £2 _ ^2 + _ A + £0 + ,e ' ^^^^ 

where we have summed over the four diagrams coming from assigning a D and a Z) to each 
initial and final line. Once again we will only need the discontinuity of the above expression 
which results in a simple expression 

Disc^(^^-^-) - i^'^i'^xi [{^i-^f-mlf' 

Next we will evaluate the XEFT loop integral in the left diagram of Fig. [71 The amplitude 
for this diagram is 

^2 r |/|2 

'^ = '^l2f^J j2^^JZ^^i^,W^'' ^^^^ 

2m^ 

where S = A — m^. The energy integral can be carried out by method of contours, the 
angular integrals are trivial, and the resulting expression is 

iM. = ID- 



--IS (-4)' 

where yU^ = 2m^(5, A is a UV cutoff and we regulate the IR divergence by inserting a factor of 
f ^/4, where f = F/ (2771,^), into the D meson propagator. In this expression we have dropped 
all terms that vanish in the A — )■ oo limit We find the imaginary part of the amplitude by 
taking the discontinuity which, using the result of Eq. ( !30ll . gives 

ImA^ = -^Disc5 -^^("a + vt^ 

3/2 



3 A4 VLl mx 
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2 ^ ( A + vri ) . (33) 



To see the effect of subtracting the above result from the result in Eq. (^^, we will consider 
the i^i — 7- A and i?2 — > A limits of the most singular terms in Eq. fl24|) (those on the last 
line of the equation), and compare it to Eq. (155]) . After converting from a decay rate to 
the imaginary part of the forward scattering amplitude, and accounting for the differences 
between the normalization of states in HH^PT and XEFT we find that the most singular 
terms in Eq. (12^ for energies near A have the form 



3 A4 nf mx Jo (p2 - ;,2)2 + ^2/4 



3 A< al mx V f 

In the first line of Eq. ( l34|) we have expressed the energy integral as a momentum integral to 
facilitate comparison with the loop integral in Eq. (l32l) . The IR divergence from integrating 
over the double pole in three-body phase space is regulated by inserting the factor 
This corresponds to exactly how the double poles were rendered finite in our calculation in 
Ref . 7| . In the second line of Eq. ( 134|) we keep only the IR divergence, the ellipsis represents 
contributions that are finite in the F — )■ limit. It is clear that XEFT reproduces the IR 
divergence in Eq. ( |24l) . The 1/F enhanced contribution from the three-body phase space 
integral cancels out of the matching calculation, and does not contribute to the coefficient 
lmA\ . Note there is a linear divergent term in the matching calculation that must removed 
by an XEFT counterterm. 

Because matching calculations are insensitive to infrared physics they can usually be 
simplified by setting any physics on order of the infrared scale to zero. In the matching 
depicted in Fig. [8] the IR energy scale is set by the typical XEFT energy scale: E 6. 
To make the IR terms in Eq. ( |2^ obvious we shift the energies: Ei = Ei — m^. Then the 
singular terms become 

tA-t ^ , (35) 

E,-A Ei-6 ^ ^ 

and the limits of integration are (i?™™ = 0,i?™'" = Qi — 2mj^). Setting the IR scale 6 to 
zero further simplifies the singular terms: 

> . (36 

Ei-6 Ei 

In addition, setting 5 — )• makes the XEFT loop integrals in Fig. |8]scaleless, and as a result 
they vanish. This greatly simplifies the matching calculation, which we do numerically. We 
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FIG. 10: A three loop XEFT diagram for T,{—Ex) with one insertion of B. 

find 

ImAf) = = -^$^-x. (105MeV)^ (37) 

X X 



where zi is given in Eq. 

The XEFT decay rate for X(3872) — )■ Xci7r°vr° can now be calculated. The dominant 

I— 5-1 ("^^ 

contribution to the decay rate comes from the diagram in Fig. Hlwith an insertion of A\ . 
The result is the same as that in of Ref. j3] but with the l/F enhanced terms omitted. In 
that paper we calculated the ratio of the branching fraction for the three-body decay to the 
branching fraction for the LO two-body decay X(3872) — )■ Xcqt^^- We now find^ 

'Br[X(3872) ^ XcivrV]\ , 

^ ^ ^ Aci J \ _ 2 9 10-3 . (38) 



Br[X(3872) XcivrO] ^ 

Because of the small propagator denominators in the graphs for X(3872) — )■ Xcitt^tt*^, this 
branching ratio is enhanced relative to the analogous branching ratio for Xco and Xc2 by 

Fig.[TO]shows a three-loop diagram contributing to H{—Ex) involving the [DDY operator 
of Fig. [9l After evaluating all three energy integrals by methods of contours, the expression 



^ A nearly identical result actually first appeared in a preprint version of Ref. [7| . The calculations of that 
version of the paper are correct even though a proper explanation of the prescription used to regulate the 
double poles is not given. 
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for the amplitude in Fig. [TO] is 



X 



X 







where / is the loop three-momentum in the left-hand loop, k is the loop three- momentum in 
the pion loop, and q is the loop three-momentum in the right hand loop. In this expression 
kr^fir^lr^qr^'j, but 171 o ^ 171^. TMs expression makes it clear that the double pole 
in XEFT is regulated by the binding momenta of the X(3872), not the width of the D*^, 
if power suppressed terms are not dropped in the propagators appearing in the loop. But 
we can also drop terms suppressed by powers of m^/m^), and the double pole will still 
be regulated by the ie prescription. Dropping mT^/rriD suppressed terms in the propagator 
denominators, we find the amplitude factors into three integrals: 

. . g'^B f dH 2fiDD* f d^q 2^dd* 



12/2m^ J (27r)^-i P + J (27r)3 q^ + 72 
X /^^2^2m^ 
J (2vr)3 " (P _ ^2 _ ,,)2 • m 

The / and q integrals yield the factor (Apos — 7)/(27r) which, when combined with the 
wavefunction renormalization, yield the terms that are identified with the nonperturbative 
matrix element in Eq. ([5]). The k integral is UV divergent and this divergence can be 
absorbed into an XEFT counterterm. The finite part of this integral is proportional to 



sj —^x^ — ie, which is imaginary. If the coefficient B is purely imaginary this cannot lead to a 
contribution to Im E(— and hence the width of X(3872). There will be a contribution to 
lmTi{—Ex) from the the real part of B. The real part of B is due to elastic D^D^ scattering 
so this contribution to lmT,{—Ex) is not related to the decay X(3872) — )■ Xcivr^vr^. It is 



ono^o 



rather a rescattering correction to the process X(3872) — t- D D 



V. CONCLUSION 

In this paper we rigorously analyzed X(3872) decays to Xcj plus pions using the OPE 
in XEFT. We reproduce our previous the results for decays to Xcj7r° for J = 0, 1, and 
2 and Xcj7i"°7r° for J = and 2. Our previous calculations Q for XciTt^tt^ suffered from 
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IR divergences due to double poles in the three-body phase space. These were previously 
regulated with the width of the D*^, T. This led terms in the decay rate that are enhanced 
by 1/r. In this paper, we argued that in XEFT these are IR divergences are regulated by 
the binding momentum of the X(3872), not the width of the D*^. In the OPE analysis of 
the present paper, we find that the IR divergences cancel in the matching calculation and 
should not appear in the expression for the decay rate. 
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